Non-Gaussian features of chaotic Hamiltonian 
transport^ 



Roberto Venegeroles and Alberto Sac 

Centro de Matematica, Computagao e Cognigao, Universidade Federal do ABC, 
09210-170, Santo Andre, SP, Brazil 

E-mail: roberto.venegeroles@ufabc.edu.br, asaa@ime.unicamp.br 

Abstract. Some non-Gaussian aspects of chaotic transport are investigated for 
a general class of two-dimensional area-preserving maps. Kurtosis, in particular, 
is calculated from the diffusion and the Burnett coefficients, which are obtained 
analytically. A characteristic time scale delimiting the onset of the Markovian regime 
for the master equation is established. Some explicit examples are discussed. 



PACS numbers: 05.45.Ac, 05.60.Cd, 05.20.-y 



I This work is dedicated to the memory of Prof. Radu Balescu. 
§ On leave of absence from UNICAMP, Campinas, SP, Brazil. 



Chaotic Hamiltonian transport 



2 



1. Introduction 

Discrete-time systems have a prominent role in many branches of nonlinear sciences. 
Hamiltonian (or area-preserving) maps, for instance, are particularly relevant for the 
modeling of classical dynamical systems [TJ [2]. The study of chaos in such systems has 
followed two main lines. The first one considers some individual trajectories in order 
to explore and characterize the system main topological properties [TJ [2]. The second 
one investigates the distribution functions of statistical ensembles or, more specifically, 
some transport properties of the associated mapsj3]. The present work belongs to this 
latter group. 

In the past, many investigations of transport properties for Hamiltonian maps 
have been motivated by the paradigmatic Chirikov- Taylor standard map|4j. Although 
considerable progress in the study of diffusion has been achieved in the last 25 years 
[3 [TJ, [91 [TOl [TTJ, H2], higher order transport coefficients have been particularly 
overlooked. The the fourth order coefficient B known as the Burnett coefficient [3] , for 
instance, plays a central role in the large deviations theory: its magnitude gives the first 
indication of the deviation of a density function from a Gaussian packet. 

Here, we consider some non-Gaussian features of the chaotic transport for 
Hamiltonian maps by means of high order corrections to the spectral properties of 
the associate Perron- Frobenius operator U. Exponential relaxation for U was rigorously 
established for hyperbolic systems by Pollicott and Ruelle [T3] IT3] . The relaxation 
rates 7 m , known as Pollicott-Ruelle resonances, are related to the poles z m = e 7m 
of the resolvent R(z) = (z — U)^ 1 . These resonances are located inside the unit 
circle in the complex z plane, despite that the spectrum of U is confined to the 
unit circle because of unitarity|15j. The normal late time evolution of density or 
correlations functions are dominated by the leading Pollicott Ruelle resonances (LPR). 
Recently, this mathematically well-established results has been confirmed in the high 
stochasticity approximation for some mixed systems [To] ITT]. In [TTJ, the LPR resonances 
are analytically calculated for the generic radial twist mapjTJ 



defined on the cylinder — n < 9 < n, — oo < / < oo. We call f{6) and a(I) = a(I + 2nr), 
respectively, the impulse function and the rotation number. The constants c, r, and K 
are assumed to be real, and K is named the stochasticity parameter. One can also 
consider non-periodic cases by taking the limit r — > oo. The LPR resonances for (TjQ) 
were obtained in [TTJ without any high stochasticity approximation, up to corrections 
of second order in the wavenumber. It is shown, in particular, that the wavenumber 
dependence of the LPR resonances determines the transport coefficients. 

In the present paper, we extend the results of [TTJ to higher order wavenumber 
corrections with the purpose of evaluating the Burnett coefficient for the map (CQ). 
Kurtosis, in particular, is then explicitly calculated. Our results are compared with 
numerical simulations for some specific models, namely the standard map, the sawtooth 



I n+1 = I n + Kf(9 n ), 

d n +l = On + Ca(I n+1 ) 



Hlod27T 



(1) 
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map, and two maps with non- linear rotation numbers: a periodic one (the tangent map), 
and a non-periodic one (the cubic map). In all cases, a very good agreement is obtained. 
Our results allow us also to infer a characteristic time scale delimiting the onset of the 
Markovian regime for the density function. We show that such a characteristic time 
scale is sharper than the others previously obtained in the literature. 



2. Statistical Analysis 

The statistical analysis of the map (0Q) is best carried out in Fourier space. The 
conditional probability that an initial state (Iq, 9q) evolves to a final state (J n , 9 n ) is 
given by / dldO p n (I, 9). The Fourier expansion of the distribution function p n (I, 9) can 
be written as 

p n (I,9) = W dqe< me +^a n (m,q) , (2) 



where the initial density is given by p — 5 {I — Io)$(0 ~ $o)> an d thus a (m,q) = 
(2ir)~ 2 ' e - l i me o+qio) _ ^he expected values of the moments I p can be calculated from the 
Fourier amplitudes a n (m,q) by 

(I?} n = (27r) 2 [(zd q ya n ( q )} q=0 , (3) 

where a n (q) = a n (0,q). An alternative way to calculate the moments is given by the 
following expression 

F n (q) = {exp[-iq(I-I )]) n , (4) 

known as the Van Hove incoherent intermediate scattering function [3]. In the limit 
q — > 0, F n (q) becomes the generating function of the generalized moments: 

Fn(q)=ex P Y, { - 1 V L C l (n) } (5) 
i=i L 

where C\ denotes the cumulant moments [1 8 [ [19]. Defining AI = I — Iq, the first 
cumulants are given by 



Ci = (AJ> n , (6) 
C 2 =((Alf)-Cl (7) 



C 3 = ((A/) 3 ) n - 3dC 2 - Cl (8) 



C 4 = ((A/) 4 ) n - 3C 2 2 - 4C 1 C 3 - QCfC 2 - C\. (9) 

The existence of cumulant moments satisfying the equation (jSj) is not assumed a priori. 
Note, however, that F n (q) is analytic around the origin q = if and only if all moments 
/(AI)^ exist and are finite. This condition breaks down in cases of distributions with 
"fat tails", like the non-Gaussian Levy stable distributions[20]. The Van Hove function 
01]) can be calculated explicitly from the density (j2J) 

F n (q) = J dldde-W-^PnW) = (2n) 2 e^ Io a n (q). (10) 
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Assuming that the evolution law for the relevant Fourier amplitude a n (q) is exponential 
for long times, 

a n (q) = exp[nr/(q)]a (q), (11) 

and that the initial relevant amplitude is given by a (q) = (2n)~ 2 e~ iqI ° , the dispersion 
rate 7(g) can be obtained from the limit n — > 00 of the the Van Hove function 

7(g) = lim -ln[F n (g)]. (12) 

Combining the generating function (jSj) with the dispersion rate (Tl2|) , one can define the 
generalized transport coefficients T> 2 i by 

1 C 2l (n) (-1)' 



V 



21 



lim — - 

71-+00 n (2l\ 



(2/)! 



dtl{q)\ q =». 



(13) 



If G\ = 0, the Einstein formula for the diffusion coefficient D is obtained by setting 
/ = 1, 



D = lim — ((A/) 5 

n^oo 2"n, \ / n 



^ 2 7 (g)| 



g=0; 



while the Burnett coefficient B is obtained for 1 = 2, 



B = lim — — 

n— >oo 4!fi 



((A/)' 



3((A/) : 



4#7(?)l 



9=0- 



(14) 



(15) 



The diffusion D and the Burnett 5 coefficients are the bases for our analysis on non- 
Gaussian features of the chaotic transport for the map ([T|). 



2.1. Accelerator modes 

Before starting our analysis, however, one should warn about the so-called accelerator 
modespQ, corresponding to fixed points (9*, J*) of (CQ): 

Kf(9*) = 2-irrL!, ca(Q = 2ttL 6j (16) 

where Lg and Lj are integers satisfying the stability condition 

\2 + cKf'(6*)a'(h)\ < 2. (17) 

Typically, trajectories diffuse normally, although some of them may be dragged along the 
accelerated modes, if they do exist. These rare events become meaningful for sufficiently 
high time scales, resulting in anomalous diffusion of Levy type for map parameters 
satisfying (IT7|) . In such a case, the diffusion coefficient behaves locally between normal 
dynamics and accelerator modes, for which one has D ~ n^ -1 for 1 < (3 < 2 [20l [2T] . In 
the case of the standard map, these divergences result in peaks for the value of D for 
K = 2rrni, with decreasing amplitude as K increases[T]. 
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3. Higher order Pollicott-Ruelle Resonances 

The dispersion rate (112D for the system (CEJ) was considered in [17 \ up to order 0(q 2 ) 
by means of the decomposition of the resolvent R(z) = (z — Z7) -1 , based on the 
projection operator techniques utilized by Hasegawa and Saphir[22] and Balescu[20] for 
the standard map. The operator U defines the law of evolution of the Fourier amplitudes, 
a n(<?) = U n a (q). Its iteration U n can be formally obtained through the identity 
§ c dzR(z)z n = 27riU n , where the contour of integration lies outside the unit circle. One 
then introduces the mutually orthogonal projection operators P = \q,0) (q,0\, which 
selects the relevant state, and its complement Q = 1 — P, leading to 

If z n 

an(q) = 77-9 dz ^5 — -»t, 1 \ a °(g)' ( 18 ) 

2m Jc z- }2j=o z ^jiq) 

where ^j{q) are the so-called memory functions for the system (pQ). The resulting 
integral is solved by the method of residues by truncating the infinite denominator 
series at j = N and then taking the limit N —>■ 00. The nontrivial leading pole was 
evaluated in the limit n — > 00 by the well known Newton-Raphson iterative method 
starting with z = 1 [17] . 

The 0(q A ) correction of the LPR resonances can be obtained by introducing into 
the denominator of the equation (fT8l) the 0(q A ) corrections to the value of z, 

z = l-Dq 2 + 0(q A ) (19) 

and repeating the same steps done in [T7]. Taking into account that ^o(q) = 1 + 0(q 2 ) 
and tyj>i(q) = 0(q 2 ), the higher order LPR resonance can be rewritten as 

00 

7 (g) = ln£(l +jDtf)* j {q) + 0(q 6 ). (20) 
3=0 

The memory functions ^j(q) are the same ones obtained originally in [TTJ for the system 
© 

Mq) = M~Kq), (21a) 
= E J-m{-Kq)J m {-Kq) G (r, mc) , (216) 



*i>a(?) = E E J- mi (-Kq)Jm ] (-Kq)Gx 1 (r,m 1 c) x 

{m} {A}t 



J 

x II QXi{ r ^ m i c ) Jrm-i- 
i=2 



i-1 



k=l 



(21c) 



The Fourier decompositions of the a (J) and /(6 1 ) functions are, respectively, 

0,(r, x) = — [d6 exp{-i[xa{r6) - W}} , (22) 

27T J 

J" m (x) = ^Jd6 exp{-t[m6 - xf(6)}} . (23) 
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Hereafter, the following convention is adopted: wavenumbers denoted by Roman indices 
can only take non-zero integer values, whereas wavenumbers denoted by Greek indices 
can take all integer values, including zero. For fixed j, the sets of wavenumbers are 
defined by {m} = {m l5 . . . ,mj} and {A}* = {Ai, . . . , Xj}, where the dagger denotes the 
restriction Yh=i \ = 0. We introduce also the following series expansion for the function 



Jm(x) = 6 mfl + °m,n X* , (24) 
n=l 

c m , n = —-. \ d6f n (6)e- me . (25) 
In n\ J 



where 



4. Transport Coefficients and Kurtosis 

The general exact diffusion coefficient can be obtained by using the definition (j!4j) and 
the LPR resonance (1201) . taking into account (12 lap - (12 left . 

D 00 

jr- = 1 + 2 a m,m Re[£ (r, mc)] a mi, mj G\i [r, m x c) x 

U 1 l m>l j=2 {m} {A}t 

i _ / k l ~ x \ 

X II ( r ' m i C ) Jm^-rm J2 X k) > ( 26 ) 

i=2 \ r fc=l / 

where the condition c 0j i = (and thus G\ = 0) is requeried, D qi = —cq^K 2 is 
the quasilinear diffusion coefficient and cr TO>m / = c^ m ^c m i cq^, [TTj . The diffusion 
coefficient (1261) gives the lowest order macroscopic description of the diffusion process. If 
the evolution process is asymptotically truly diffusive, then the angle-averaged density 
should have a Gaussian contour after a sufficiently long time. A first indication of the 
deviation of a density function from a Gaussian packet is given by the Burnett coefficient 
B defined by ( Tl~5l) . The dimensionless fourth order cumulant 



«(„) = <^>f (27) 

{(Ami 

is usually called the kurtosis. For Gaussian densities, B = and k = 3 for all times. 
Combining (fI4|) . (fl5l) . and (f27|) . we obtain for sufficiently long times 

K = 3 + ^!. (28) 

The Burnett coefficient I? can be evaluated in the chaotic regime by truncating of the 
resonance (1201) at j = 2 so that (fl5l) can now be applied to (1201) . yielding 

5 -D 2 + 2DD ql - 1 -Dd^ 1 + ± 

calculated at q — 0. In order to verify these results, we calculate the Burnett coefficients 
and its respective kurtosis for some particular cases of maps (CQ) and compare with the 
respective numerical simulations. To this purpose, it is important to choose intermediate 
values of n, since exaggerated values tends to wash non-Gaussian fluctuations away 
quickly. 



B « --D 2 + 2DD gl - -D d 2 ^ 1 + -d 4 q ^ + + tt 2 ) (29) 
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4-1. The standard map 

The standard map corresponds to the choice ca(I) = I and f(8) = sin 8 in (OQ). The 
memory functions (12 1 a[) - ( f2TcT) are rather simple for the case of linear rotation number 
since one has Q\(l, x) = 5\ jX . For the standard map, J mix) is the Bessel function of the 
first kind J m (x), D q i = K 2 /4 and o- m ,m' = (±#m,±i)(±£m',±i)- The final expressions for 
D and B in the chaotic regime for the standard map are given by 

^-=l-2J 3 {K)+0{ir 1 ), (30) 
= - \ + UK) + 2J 2 (K) + UK) + ^U2K) + OiK- 1 ). (31) 
The kurtosis ( |28i) can be evaluated straightforwardly. Fig. [1] depicts the calculated 
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Figure 1. Kurtosis as a function of K for the standard map for c = 1 and 11 = 10. The 
dots correspond to numerically computed values and the solid line to the theoretical 
result in the chaotic regime, up to terms of order OiK^ 1 ). For each value of K, 10 5 
random initial conditions are considered in the numerical simulation. This plot shows 
good agreement with numerical calculations although the accelerator modes give rise 
to spikes in the figure. The dashed line corresponds to the Kurtosis calculated from 
the results of Tabet et. al presented in [21]. See Sect ion [431 for further details. As one 
can see, the present results give rise to a superior description. 

kurtosis and the results of numerical simulations for the standard map. We notice, in 
particular, the presence of accelerator modes. 
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4-2. The sawtooth map 

The sawtooth map corresponds to ca(I) = I and f(6) = 9. As for the standard map, 
^(1)^) — $\,x- On the other hand, in this case, J m {x) = sm }j^x)^ , D q i = K 2 tc 2 /6, and 

s? / -i\m, — m' 

o- m , m > = rn m , , leading to 

= 1 - WlC), (32) 

D q l 6 

= - 5 - + YTY^ K) + (2 + i^ S *™, (33) 

both up to (D(K~ 2 ) order, where the function Sj(K) and Cj(K) are given by 

- 72sin(7rmif) 

^2j+l(A) = 2^ fir , U^2j+l ' ( 34 ) 



m= 

oo 



i (K + 2)(7rm) 2 



r , _ ^ 72 cos(TrmK) 

and S2j(K) = C2j+i(K) = 0. Fig. [2] presents the comparison of the kurtosis calculated 
from fl32l) and (1331) with the numerical simulations for the sawtooth map. As one can 




Figure 2. Kurtosis as a function of K for the sawtooth map for c = 1 and n = 10. 
For each value of -fT, 10 5 random initial conditions are considered in the numerical 
simulation. This plot shows excellent agreement between the theory and the numerical 
simulations. 

see, the absence of accelerator modes contributes with the excellent agreement of the 
numerical simulation with the theoretical predictions. 

A pertinent comment here is that, as one can see from Figures [1] and EJ for a fixed n, 
the limit for the kurtosis k is not 3 as K — > oo. From the asymptotic values of ( l30l) - (!3Tl) 
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and (13"2"j) - (l3"3"j) . it is easy to show that, for n = 10, the kurtosis in the limit K — > oo 
tend to the values 57/20 and 72/25, respectively, for the standard and sawtooth maps. 
Similar results hold also for other maps. The conclusion that, for fixed n, the limit of 
high stochasticity is not enough to assure a Gaussian regime is interesting and certainly 
would deserve a deeper analysis. 



4-3. The tangent map 

We call the tangent map the choice of f{6) = 9 and a(I) = tan(J/2) in ([1]). It is our 
first example of a map with a non-linear rotation number. The functions J m {x) are the 
same ones from the sawtooth map. The functions Qe(l, x) can be calculated by complex 
residues from their definition (122j) . Introducing the variable s = tan(#/2) and taking 
into account the identity iarctans = arctanhzs, we have 

For £ = 0, ( |36i) has two single poles on the complex plane located at z = ±z. For positive 
x and negative x, one closes the integration path of (I36I) . respectively, by the positive 
^s{z) and negative $s(z) semiplanes, giving simply 

o (l,aO = e-H. (37) 

Let us consider now £ > 0. Notice that ^(1, x) can be evaluated for negative £ by 
observing that Q-e(l, x) = Qe(l, —x). The integral fl36l) has a unique pole of order £ + 1 
at z = i for £ > 0. In such a case, for negative x we can close the integration path in 
the negative $t(z) semiplane and conclude that Qg(l,x) = for £ > and x < 0. For 
positive x we close the integration path in the positive ^s(z) semiplane and obtain 

_f-11 



Gi(l,x) 



+ . (38) 

J s=l 

for £ > and x > 0. As illustrative examples, the very first functions fl38l) are given by 



0i(l,x) = 2a;e- x , (39) 

&(l,a;) = 2x(x- l)e~ x , (40) 

3 (l,x) = ^x(2x 2 -6x + 3)e" x , (41) 

for x > 0. 

The expressions for D and B are obtained from (|26l) and (1291) . respectively. Fig. 

|3] depicts the comparison between the calculated diffusion and kurtosis for the tangent 
map and the numerical simulations. 



4-4- The cubic map 

We call the cubic map the choice oc(I) = J 3 /3 and f{6) = sin6 in ([1]). It is our example 
of map with non-periodic non-linear rotation number. Due to the non-periodicity, some 
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Figure 3. Diffusion (a) and Kurtosis (b) as a function of K for the tangent map with 
c = 3 and n = 10. The dots corresponds to the numerically computed values and 
the solid line to the theoretical result in the chaotic regime. For each value of K, 10 5 
random initial conditions are considered in the numerical simulations. Both plots also 
exhibit good agreement between the predictions and the numerical simulations. 

parts of our formalism shall be modified. First, we rewrite the function (122j) in the 
following way 

g x (r,x) = r due- i(xa(u) - Xr ' lu) . (42) 

27ir J-nr 
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The non-periodic rotation number corresponds to the limit r — > oo. In this case, we 
apply the following overall replacements in our formulas 

r- l \ — ► s, (43) 

^2 — > / dsi x . . . x dsj , (44) 



r 



-j 



{X 1 ,...,X J } 

r Q\{r, x) — > G(s, x), (45) 

where the function Q(s,x) is now given by 

g(s,x) = — [ dl e -<Mi)si) . ( 46 ) 

Performing analogous calculations to the periodic case, we obtain the new set of memory 
functions 

Ml) = M-Kq) , (47a) 

=E^-™ (-Kq)J m (-Kq)g(0,mc) 7 (476) 



*i>2(g) = ^2J-m 1 (-Kq)Jm j (-Kq) J d8Q{sx,mxc) 



{m} 

j 



X 



J\g{s h m i c) 



i=2 



X 

St 

/ »-l 



k=l 



(47c) 



where cZs" = ds\ x . . . x ds_j, with the domain of integration defined by = 
|(si, . . . , Sj) : Y^{=i s i — 0|. Finally, the diffusion coefficient will be given by 

j-y oo oo „ 

— = 1 + 2 ^ cr m , m Re[£?(0, mc)] +51X1 ^.mj / dsG(s l} m x c) x 

-^V m=l i=2{m} 1/51 

i=2 \ fc=l / 

(48) 

Note that, for non-periodic linear rotation numbers, we have g(s,x) = 5(s — x), and 
the diffusion formula (148]) coincides with the periodic linear rotation number version of 

(EE}. 

For the cubic map, the function £?(s, 2) can be calculated by means of Airy 
functions [23] 

g(s, x) = x- 1/3 Airy(-x" 1/3 s). (49) 

The power dependence x~ 1//3 in ( 149]) may create a false impression of divergence of the 
series (148]) for the cases where c 1//3 <C 1. In order to avoid this problem, one can define 

c = cic 2 , c^ 3 I = x, (50) 

and rewrite the cubic map as 

x n+ i =x n + Kc 2 /3 sen 9 r , 



'n 1 



(51) 

6 n+1 = 6 n + cix 3 +1 /3 mod27r. 
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From (jSDJ), it follows that (Ax) 2 = c 2 /3 (AI) 2 and D { f 
D/D q i for the cubic map is invariant under the rescaling 



2/3 n (/) 



D 



TT ( c > *0 = TT ( c i c 2 5 *Q = 7T ( c i> Kc : 



1/3 > 



9/ 



41 



Assuming that 



K > c 



-1/3 



the rescaling 



ci = if 3 c, 



C 2 = 7T 



12 

Hence, the rate 
(52) 

(53) 
(54) 



prevents any potential problem of divergence for maps with small c. 

For small values of K, the memory functions ^j>2 (<?) give rise to high oscillatory 
combinations of Airy and Bessel functions, with integrals that are very difficult 
to estimate. On the other hand, the high stochastic condition (1530 implies that 

^ In 

Airy( Cl 1 s) ss Airy(O) even for sufficiently high values of \s\ for which Bessel functions 



already decay as J m {s) ~ e 
map can be estimated by 

D 



-1/2 



D, 



1 



ql 



<E 



{ (Re 1 / 3 )!' 



. Thus, in this regime, the rate D/D q i for the cubic 



(55) 



where 



b 



2Airy(0), 



(56) 



J>2 



W(0) ( ? ( Kx 



(±5 mi ,±l)(±5 m ,,±l) 



x 



x / dsi x 



■[2Aiiy(0)]' £ 



(±<5, 



mi,±l 




(777.1 x . . . x rrij) 



Airy 2 (0)<5 ii2 , (57) 



{mi=2Zi+l} 

leading, as expected, to Wxuk^ooD = D q i. Indeed, the quasilinear regime for the 
diffusion is rapidly attained for high values of K, without oscillations, see Fig. HI 
Similar results hold for any rotation number of the type cx(I) oc I p , for p > 1. In such 
cases, Q(x, 0) oc and ( |53j) shall be replaced by K ^> c~ x l' p . 

Regarding the kurtosis for the cubic map, similar arguments can be used to show 

that 

k ~ 3 - — , 

on 

also without oscillations for K ^> c -1 / 3 . Since the diffusion is quasiliner for such values 
of K, the memory functions *f?j can be disregarded for all j > 1. In this limit, from (|2"9|) . 
we have simply B/D 2 ~ —1/5, and, hence, leading to fl58|) . Fig. [5] depicts the behavior 
of the kurtosis for the cubic map. 



(5f 
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Figure 4. Numerically calculated diffusion for the cubic map with c = 10~ 4 (a) 
and c = 1CP 1 (b). For each value of K, 10 5 random initial conditions are considered 
in the numerical simulations. According to our results (see Section |4.4[) , quasilinear 
behavior is expected to occur for K ^> c -1 / 3 . From the figures, quasilinear diffusion 
clearly takes place for K > 21 and K > 2, respectively, in accordance to our theoretical 
predictions. 



4-5. Comparison with previous approaches 

It is instructive to compare our results with others previously obtained in the literature 
in the context of standard map. Tabet et et/[2l], for instance, do not calculated the 
Burnett coefficient explicitly. They used, instead, the Fourier path technique [U [5] to 
calculate the moments ((A/) 4 ) n and ((AJ) 2 ) n , from which we can evaluate the respective 
rate B/D 2 by means of the definitions (114p and ( fi"5{) . One has 

^2 = "7 + 3J 2 ( R ) + 2J 4 (2K) + 0{K- 1 ). (59) 

On the other hand, Balescu[25j utilizes a related, but somewhat different approach from 
ours to calculate kurtosis for the standard map. He also does not calculate explicitly 
the kurtosis nor the Burnett coefficient, but he present a non-Markovian approach for 
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Figure 5. Numerically calculated kurtosis for the cubic map with c = 10 (a) and 
c = 10 _1 (b). For each value of K, 10 5 random initial conditions are considered in 
the numerical simulations. As the diffusion coefficient (Fig. |4j, quaselinear behavior 
is expected to occur for K > 21 and K > 2, respectively. 



the relevant amplitudes, for which the propagator F n (q) is given by 

F n « % + (n - 2)^- 3 ^ 2 + l(n - 5)(n - 6)*J" 6 *1 ( 60 ) 

for ri > 6, where memory functions ^j>s contributing with terms of order 0(K~ l ) are 
disregarded. The corresponding Burnett coefficient can be calculated by means of (fT2!) 
and (IT5]) . leading exactly to the same result we have gotten here 

Jj2=~\ + MK) + MK) + UK) + i J 4 (2K) + O^ 1 ), (61) 

compare with (1HT1) and (I5UI) . Although the equations (1551) and (IBTl) have the same 
asymptotic value, they are enough different even in the high stochastic regime. However, 
it becomes difficult to note differences in the kurtosis calculated in both cases for large 
values of n (Tabet et al used n = 50, for instance). For relatively smaller values, the 
equation (16T1) gives a much better agreement with numerical simulations than (l59l) . as 
one can see in the FigJTJ 



5. Gaussian characteristic time scale 



Equation (1251) suggests the existence of a characteristic time: the Gaussian time scale 
no, defined by 



Tin, = max 

K 



(62) 
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where \x\ = inf{n e Z\x < n} is the ceiling function. For n 3> n G , the transport 
process is typically Gaussian, up to order 0(q 6 ). In fact, for n ^> n G , the expansion of 
the propagator exp [77,7(5)] gives the well known Gaussian density 

exp[n 7 (g)] = 1 - Dnq 2 + \Dn 2 ^- g 4 + 0(q e ) 

« exp(~Dnq 2 ) + C(g 6 ), (63) 

where k is given by equation fl28l) . The Gaussian time scale n G has also a second and 
no lesser important interpretation: in the regime n > ng, the time evolution of the 
relevant amplitudes becomes Markovian. This is easy to realize since, for n ^> n G , the 
higher order corrections can be neglected and the LPR resonance ( )20|) . which gives the 
exact expression for the normal diffusion coefficient [IT] , can be taken as the approximate 
propagator. Hence, for n ^> n G , equation ffTTl) holds perfectly and, furthermore, becomes 
the following Markovian master equation 

*»(<?)« ( 64 ) 

The Gaussian time scale n G is sharper than some characteristic times obtained 
previously in the literature, as, for instance, Balescu's memory time hm introduced 
in [25]. Balescu obtained the propagator (I60I) by means of the general non-Markovian 
Bandtlow and Coveney master equation[26j. The convolution of the master equation is 
truncated at the memory time n M , leading to 

n c 

a n +i (q) ~ ®j (?) a n-j (?) , (65) 
3=0 

where 

j n for n < n M , 

n c = \ (66) 
I Um for n > n M - 



Notice that, for n 3> n M , the equation (j65]) becomes the Markovian equation 
By means of some numerical experiments with the decay of memory functions for the 
standard map, Balescu concludes that tim = 4 and, in such a case, obtained equation 

dSOD. 

Evidently, we expect that the two time scales to be related by n G — Um + 1- The 
advantage of n G is that it can be calculated judiciously by means of fl62l for a generic 
class of systems like (CE]). Moreover, typically, n G is sharper than Balescu's memory 
time. In Fig. [6] we can observe the behavior of the rate [6|5|/D 2 ] as a function of 
K for the standard map and for the sawtooth map. The theoretical prediction given 
by (1621) is confirmed for K > 4 for the standard map (K > 2 for the sawtooth map). 
However, it is important also to point out that the Gaussian time (1621) evidently fails in 
the description of the weak-stochasticity regime, since the effects of the KAM surfaces 
and the stable islands become increasingly important as K decreases, requiring not 
only the calculation of further memory functions but the inclusion of the source term 
involving initial conditions (see [26]). For instance, let us consider the equation (l65il for 



Chaotic Hamiltonian transport 



16 



CNI 

Q 

CD 

c 

q3 
o 



Q 

m 

c 

O 



6j 
5 

4- 
3- 
2- 
1 - 



(a) 





20 40 60 80 100 120 




K 



Figure 6. Gaussian time scale as a function of K for the standard map (a) and for the 
sawtooth map (b), both with c = 1 and n = 10. For the standard map, we have fig = 3 
(or riM — 2) for K > 4 (for the sawtooth map, for K > 2), in excellent agreement with 
our expectations. For larges K, n^j falls to 1 indicating the absence of correlations of 
order K~ x l 2 for large n, for both cases. 



the general linear rotation number case by applying the Balescu approach for um — 2 
(see Appendix). New calculations of the Burnett coefficient by means of the propagator 
( 1A.6I) gives 



£ ~ + 2£>D 9J + -c\ 4 (v|> + H/ 2 ) 9=0 , (67) 

that is exactly the equation (1291) for = 0! Hence, we should expect in the chaotic 
regime 

n G = 3 for ca(I)=I. (68) 



6. Summary and Discussion 



In this paper, we have performed a detailed analysis of the non-Gaussian aspects of the 
normal transport in Hamiltonian discrete systems. The general class of area-preserving 
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maps represented by (pQ) has been chosen because it comprises a large number of physical 
situations and has the paradigmatic standard map as a particular case. The map (JTJ) 
was recently studied in [TTj , where the basis for the study of the higher order transport 
coefficients, including nonhyperbolic systems, was initiated. The LPR resonances of the 
system (CD) were enhanced in its wavenumber dependence with corrections of order g 4 , so 
that the corresponding fourth order Burnett coefficient could be evaluated. Numerical 
simulations were performed for four particular cases of (TjQ) and excellent agreement with 
the theoretical predictions is obtained. 

We have established also a Gaussian time scale no given by equation (|62|) . For 
n ^> ric, the density function acquires a Gaussian countour and its time evolution is 
ruled by a Markovian master equation. We also show that uq is related with the memory 
time riM defined by Balescu [25] by nc = Um + 1- For maps such that ca(I) = I we 
conclude that tlq = 3, a sharper result than Balescu's memory time tim- 
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Appendix A. Balescu's approach for tim = 2 

In this appendix we show that the cuttof time tim in the Balescu equation ( l65i) can 
be taken as tim = 2 in the chaotic regime. For every calculation, we retain arbitrary 
powers of \I/o and terms of order ^2 and (recalling that ^1 = for the linear rotation 
number case). Equation ( 1651) in the form a n+ i(q) = F n (q) a n (q) where 

n c 

F n+ i(q) = E *i(?)^n-i(g), ^o(g) = 1, (A.l) 

3=0 

with n c = n for n < Um and n c = Um for n > tim- The next six propagators F n are 
given by 





= ^0 




F 2 


= < 


F 3 




+ ^2, 


F 4 




+ 2^0*2, 


F, 


= *o 




F 6 




+ 4*g* 2 + 



Thus, we can write a general expression for F n : 

F n = % + ar n *r 3 *2 + J/n*?~ 4 *l, (A.2) 
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with the following initial conditions (note that ^0^2 ~ ^2 ~ ? 4 f° r finite k) 

\ Vn<5 = 0, 2/6 = 1- 

On the other hand, the equation (IA.10 gives 

F n = #0^-1 + ^2^-3, n > 3. (A.4) 
Substituting the equation ( 1A.2I) into ( 1A.4I) and comparing the coefficients we obtain: 

{•^n %n— 1 "I - I? p.^ 



Solving the system flA. 5|) with the initial conditions (1A.3|) we finally obtain the 



propagator F n for n > 6: 



F n = *J + (n - 2)^~ 3 ^ 2 + i(n - 5)(n - 4)^" 4 ^, (A.6) 
which may be compared to fl60l) . 
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